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ABSTRACT 



Context. We analyse the line and continuum spectra of the symbiotic system CH Cygni. 

Aims. To show that the colliding-wind model is valid to explain this symbiotic star at different phases. 

Methods. Peculiar observed features such as flickering, radio variation, X-ray emission, as well as the distribution of the nebulae and 

shells throughout the system are investigated by modelling the spectra at different epochs. The models account consistently for shock 

and photoionization and are constrained by absolute fluxes. 

Results. We find that the reverse shock between the stars leads to the broad lines observed during the active phases, as well as to radio 

and hard X-ray emission, while the expanding shock is invoked to explain the data during the transition phases. 

Key words, binaries: symbiotic - stars: individual: CH Cyg 



1 1. Introduction 

' The symbiotic star (SS) CH Cygni (CH Cyg) has been observed 

since the late '80s of the XIX century. It was identified as an 

anonymous red semiregular variable by Gaposchkin's 1952 anal- 

. ysis of Harvard patrol plates, and classified as a M6-M7 spec- 

' trum star with a 90-100 day period and about one-magnitude 

■ amplitude. This classification seemed to be suitable until 1963 
when a hot blue continuum appeared, along with a late-type 
spectrum and emission lines of H, He, [Fell] and [Call], evi- 

, dence of a composite nature. This drove Deutsch (1964) to rec- 

■ ognize in CH Cyg the brightest SS at visual wavelengths. Since 
' then, active and quiescent phases have been alternating without 

a predictable trend. 

■ At a distance of 270±66 pc (according to the Hipparcos par- 
allax measurements - Munari et al. 1997) CH Cyg is one of the 
most studied, yet poorly understood SSs. It has been observed in 
the last thirty years both photometrically and spectroscopically 
from radio to X-ray wavelengths (Kenyon 2001). 

The basic configuration of the CH Cyg system is still de- 
bated, despite the large amount of multifrequency observations 
carried out since the first recorded symbiotic activity. In partic- 
ular, it is discussed whether it is a binary or a triple star, and, in 
the latter case, whether composed of a red giant/white dwarf pair 
with a G star companion (Hinkle et al. 1993) or a red giant/white 
dwarf pair with an additional red giant companion (Skopal et al. 
1996). The explanation of the total hot component luminosity in 
a binary star scenario might be a crucial problem: specifically, 
the suggested binary orbital parameters (e.g. Mikolajewski et al. 
1992) would imply an accretion rate onto the white dwarf not in 
agreement with that requested by the observed hot component 
luminosity. On the other hand, adopting the binary star scenario 
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the orbital solution of the triple model (Hinkle et al. 1993) would 
require a low inclination angle. This is in contrast with the jet ap- 
pearing in the plane of the sky, implying an accretion disk and/or 
an orbital plane virtually edge-on. Furthermore, the low inclina- 
tion implied by Hinkle et al. (1993) would result in an unreal- 
istically low mass (0.2 M© ) white dwarf. Finally, the similarity 
of the 756 day period of CH Cyg with non-radial pulsations in 
a few semi-regular variables drove Schmidt et al. (2006) to in- 
terpret this periodic trend more as a non-radial pulsation of the 
cool component rather than a real orbital period. 

Several theoretical papers have appeared presenting SSs 
as colliding-wind (CW) binary systems (e.g. Kenny & Taylor 
2005). The observational evidence that actually both the com- 
ponent stars of a SS lose mass through strong stellar winds 
(Nussbaumer et al. 1995), along with the presence in CH Cyg 
of radio, optical and X-ray jets (Karovska et al. 2007), have con- 
firmed the importance of shocks for a real understanding of such 
outstanding system. 

CH Cyg is characterized by the presence of an eruptive white 
dwarf, showing episodes of accretion fed by the high-density 
wind coming from the M star, and of wind-wind interaction lead- 
ing to shock waves. The nebulae downstream of the shock fronts 
show the spectroscopic marks of a turbulent kinematics, result- 
ing in an inhomogeneous and highly disrupted envelope. 

In this paper we review the recent history of CH Cyg, in 
the light of a colliding wind scenario (e.g. Girard & Willson 
1987, Kenny & Taylor 2005) that has been used to explain many 
SSs (e.g. Formiggini, Contini, & Leibowitz 1995, Contini 1997, 
Contini & Formiggini 1999) and recently updated including a 
detailed analysis of dust features (Angeloni et al. 2007a,b,c). 

We interpret the observational data throughout both the ac- 
tive and quiescent phases, from the 1978 great outburst up to 
date, with the aim of investigating the physical and morphologi- 
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Fig. 1. AAVSO light curve (10-day means of observations) of 
CH Cyg from 1948 to 2000. 



cal conditions of the system by modelling the line and continuum 
spectra emitted from the different nebulae. 

We adopt the binary scenario, considering the ~ 756 day pe- 
riod as due to the giant pulsation mode, and the ~ 15 year period 
as the actual orbital period. 

Following the method used in previous SS investigations, we 
analyse the continuum spectral energy distribution (SED) in dif- 
ferent epochs, instead of the light curve through the years. At 
each epoch we model the data collected from radio to X-ray 
constraining the results by the analysis of the line spectra. The 
continuum accounts for bremsstrahlung and dust reprocessed ra- 
diation from the shocked nebulae, radiation from dust shells, as 
well as radiation from the stars. Our method is particularly suit- 
able for disentangling the multiple contributions of gas and dust 
in different conditions at different phases. 

We focus on the nebulae within the shocked wind-wind in- 
teraction regions and on the shells ejected by the giant star, con- 
straining the models by the observed absolute fluxes. 

A short summary of observational history of CH Cyg is pre- 
sented in Sect. 2. The colliding wind model and our application 
to CH Cyg are presented in Sect. 3. The next sections are devoted 
to the analysis of the individual active and quiescent phases. As 
CH Cyg has not been observed regularly in the different wave- 
length ranges, except of the optical one, the SED profiles in some 
epochs present some data gaps: for instance, there are no UV 
and X-ray observations for the 1988-1990 quiescent phase, nor 
mid-IR data for the 1998-2000 active phase. This led us to focus 
on specific physical processes in each phase: the basic models 
which result by disentangling the different contributions to the 
SEDs and to the line ratios at different days during the active 
phase 1978-87 are presented in Sect. 4; in Sect. 5, different types 
of IR obscurations occurred during the 1988-90 quiescent phase 
are indicated; near-UV and IR variability during the 1991-95 ac- 
tive phase are discussed in Sect. 6, while Sect. 7 (1996-97 quies- 
cent phase) deals with the complex nature of the radio emission. 
The very expanded nebulae and the X-ray emission are modelled 
throughout the 1998-2000 latest active phase (Sect.8) up to date 
(Sect. 9). Concluding remarks and a summary of the results are 
presented in Sect. 10. 



2. Observational progress 

The symbiotic classification was introduced on the basis of op- 
tical spectral features, however, a real understanding of SSs, and 



CH Cyg in particular, has been possible only by space-based as- 
tronomy. 

A key-role was played by the long-lived lUE, which opened 
the UV window allowing to monitor the temporal evolution of 
many emission lines arising from the nebulae and the continuum 
coming from the hot star (see Sect. 4 for a detailed modelling of 
such temporal evolution). This great effort revealed that CH Cyg 
is much more complex than thought before. Its puzzling nature 
triggered a series of new radio observations (e.g. Taylor et al. 
1985, 1986) attempting to understand the actual origin of the 
nebular emission and to constrain a self-consistent view in the 
light of the most accepted scenario at the time. It was concluded 
that the symbiotic behaviour was due to the dense cool wind 
photoionized by the hot star radiation (e.g. Taylor et al. 1988). 

The optical and near IR spectrophotometry throughout the 
years (e.g. Taranova & Yudin 1988, 1992; Munari et al. 1996; 
Skopal et al. 1996, 1997, 2000, 2002, 2004, 2007) have aUowed 
to follow small changes in the emitting properties of CH Cyg, 
without leading to a substantial breakthrough. 

Unexpected results came from the first X-ray observations 
(e.g. Leahy & Taylor 1987) that for the first time made clear 
that the energetic of the whole system was hard to explain in 
term of standard photoionization models. It became evident that 
shocks, and in general hydrodynamic phenomena (e.g. fluid in- 
stabilities), play a primary role in shaping the energy distribution 
along the whole spectral interval. Further observational evidence 
of stellar wind collision in CH Cyg derived from the most recent 
X-rays observations (Ezuka et al. 1997, Galloway & Sokoloski 
2004, etc). 

The IR spectral range is not yet fully explored. It was sug- 
gested that "symbiotic dust" carries as much information as sym- 
biotic emission lines (Angeloni et al. 2007c), and the unique 
CH Cyg mid-IR spectra recorded by the ISO satellite (Schild 
et al. 1999) complicated furthermore the system overall pic- 
ture. Multiple dust shells, mixed grain chemistry, equilibrium 
between grain formation and destruction processes in the turbu- 
lent nebulae (Angeloni et al. 2007a) are worth investigating with 
the new era of mm-IR facilities, such as ALMA and Herschel. 

Fig. 1 shows the light curve of CH Cyg throughout many 
years. The first ever recorded symbiotic activity of CH Cyg 
started in September 1963 and was over by August 1965, fol- 
lowed by an other period of increased activity in the years 1967- 
1970. Then, CH Cyg experienced a dramatic eruption in 1978 
when it became as bright as to be visible by unaided eye (visual 
mag. 5.6). This great outburst ended in 1986 with the emission 
of bipolar radio and optical jets (Solf 1987) and by a decline 
of about 2.5 visual magnitudes (Taylor et al. 1986). After 1986, 
the decline in magnitude continued, even though not regular nor 
smooth. Meanwhile, shorter and dimmer outbursts (1991 - 1995, 
1998 - 2000) have been accompanying the fading trend (Fig.lTJ. 
Between these periods of activity the spectrum has exhibited al- 
most exclusively red giant characteristics, and U B V light curves 
with significant variations at a level of about 1.5 mag. 

Since the beginning of the 2000, CH Cyg persists in a qui- 
escent phase at rather bright magnitude, despite a 2 mag decline 
in the U B V light curves during the 2006 June-December pe- 
riod (Skopal et al. 2007). An apparently anomalous extreme IR 
fading has been recently pointed out by Taranova & Shenavrin 
(2007). 

In the Appendix, we present the references to the observa- 
tional data we collected from the literature. In the quoted works 
the reader can find the journal of observations, so in the Tables 
we have indicated only the spectral domain and the year, ad- 
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Fig. 2. A very schematic picture of the CW model and the dust 
shells (adapted from Kenny & Taylor 2005). 



dressing to the specific reference for any further information or 
technical detail. 



3. Theoretical framework 

3.1. CW models for SSs: a short review 

The first models for the calculation of radio emission from SSs 
adopted a single star picture with a totally ionized wind (e.g. 
Seaquist & Gregory 1973) Then, Willson et al. (1984) repro- 
duced the line profiles in terms of the interaction of a low veloc- 
ity wind from the red giant with a high velocity wind from a WD 
in a detached binary system where both stars have spontaneous 
winds. The line forming regions are located in different portions 
of the interaction shell and in the two winds. 

A binary system composed of a cool, mass-losing star and a 
hot companion was also adopted by Taylor & Seaquist (1984), 
reproducing the radio spectrum by a simple steady-state model 
in which a portion of the stellar wind is photoionized by the hot 
companion. 

Girard & Willson (1987) presented the colliding wind model 
which is at the basis of further developments. Namely, the winds 
from a late giant and from a hot WD collide between the stars 
and material is accumulated along the boundary, forming an ini- 
tially spherical shell that is deformed by the wind of the primary 
(giant star): eventually, the shell reaches a steady state configu- 
ration. A large amount of turbulence perpendicular to the mean 
flow is expected. 

Both in the steady state formulation and in the dynamic shell 
model the shell is sufficiently thin to be considered having zero 
thickness. 

The colliding wind model of Kenny & Taylor (2005) dif- 
fers from previous ones because assumes mixing of shocked 
material from both incoming streams instead than postulating a 
self-sustaining contact discontinuity. The model extends the CW 
model by the derivation of an adiabatic temperature, the consid- 
eration of radiative cooling, the inclusion of thermal pressures 



in the incoming winds, and the treatment of interaction shells of 
finite thickness and density. They succeed in reproducing radio 
images and continuum spectra. 

Collision of the winds and the derived shock network have 
been invoked to resolve specific discrepancies between observa- 
tions and theory during the active phases of a symbiotic phe- 
nomenon. Nevertheless, in order to simultaneously reproduce 
the UV, optical, and IR line spectra as well as the absolute 
fluxes of the continuum SED from radio to X-ray (Contini 1997, 
Contini & Formiggini 2003, Angeloni et al. 2007 a,b), the shocks 
cannot be neglected even during the quiescent phases. 



3.2. A schematic modei adapted to the detailed caicuiation 
of the spectra 

Recent works on a significant sample of dusty SSs have been re- 
vealing a common trend in the IR spectrum properties, such as 
similar, recurrent shell temperatures and unexpected grain chem- 
istry (Angeloni et al 2007b, 2008 in preparation). 

In order to take into account the importance of dust in SSs, 
our model adds to the colliding wind binary system the dust 
shells emitted by pulsation of the cool star A very schematic 
picture of the system is shown in Fig. |2] In this figure the ba- 
sic shock fronts appear The dynamical effects of instabilities, 
such as fragmentation, will be explained in the text, as well as 
evaporation and sputtering of dust grains in the shells. 

The shocks associated with the colliding winds are strongest 
along the line connecting the two stars. The interface between 
the two winds contains two shocks: one facing the fast wind, 
with Vs~ 1000 km s 'but low mass flux, and one facing the slow 
wind, with Vs~ 30 km s 'and considerably higher flux corre- 
sponding to mass loss rates of ~ 10"'' or 10"^ Mq year"'. To a 
first approximation they can be considered standing shocks. We 
will refer to the strong shock facing the WD. In the expanding 
region around the hot star away from the direction towards the 
cool star, we will refer to the expanding shock outward the sys- 
tem. 

We are interested in the calculations of the line and con- 
tinuum spectra throughout the different epochs, so we refer to 
schematic nebulae: downstream of the strong standing shock be- 
tween the star, downstream of the expanding shock, and down- 
stream of the dust shell edges. The radiation emitted from the 
different nebulae and dust shells are summed up with different 
weights which depend on the orbital phase, the pulsation period 
of the red giant, and on activity of the WD. 

The nebulae downstream of the shock fronts are compressed 
by the shock and heated by both the shock and the photoion- 
izing flux from the WD. Downstream of the reverse shock, the 
radiation flux from the WD reaches the shock front edge of the 
nebula, characterized by very high densities and by a high ion- 
ization parameter, being the shock front close to the WD. The 
nebula downstream of the shock expanding outwards the sys- 
tem is characterized by lower densities. In this configuration, the 
photoionizing flux from the WD, with an ionization parameter 
diluted by distance, reaches the edge of the nebula opposite to 
the shock front. Furthermore, the photoionizing flux may be ob- 
structed by some intervening matter before it reaches the outer 
nebula, hence leading to shock-dominated spectra, i.e. showing 
only the characteristics of the shocks. 
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3.3. The SUM A code 

The line and continuum spectra (free-free and free- 
bound radiation, hereafter bremsstrahlung) emitted by 
the gas downstream are calculated by SUMA (see 
http://wise-obs.tau.ac.il/~marcel/suma/index.htm for a de- 
tailed description), a code that simulates the physical conditions 
of an emitting gaseous cloud under the coupled effect of 
photoionisation from an external radiation source and shocks, 
and in which both line and continuum emission from gas is 
calculated consistently with dust reprocessed radiation in a 
plane-parallel geometry. In particular, dust grain heating and 
sputtering processes are also included, as well as free-free 
self-absorption in the radio range, so evident in SSs (e.g. Hl-36, 
Angeloni et al. 2007b). 

Before proceeding with the modelling of the spectra, we 
would like to point out that the plane-parallel condition is an ap- 
proximation which is a necessary constraint in composite mod- 
els (shock -Hphotoionization). However, the plane-parallel condi- 
tion is roughly valid in the interbinary region where the collision 
of the winds is head-on. It is valid in the circumbinary region 
where the radius of the expanding shock created by the head-on- 
back interaction of the outflowing winds from the stars is large 
enough. We have adopted it also for the dusty shells because they 
can approach large radii. 

The input parameters are: the shock velocity V,, the 
preshock density no, the preshock magnetic field Bq, the colour 
temperature of the hot star T,, the ionization parameter U, the 
geometrical thickness of the emitting nebula D, the dust-to-gas 
ratio d/g, and the abundances of He, C,N, O, Ne, Mg, Si, S, A, 
Fe relative to H. Solar abundances are adopted (Allen 1973). 
The distribution of the grain radius downstream is determined 
by sputtering, beginning with an initial radius of 0.2 //m. 

The calculations start at the shock front where the gas is 
compressed and thermalized adiabatically, reaching the max- 
imum temperature in the immediate post-shock region. The 
downstream region is cut in many plane-parallel slabs (up to 
300) with different geometrical widths calculated automatically 
in order to account for the temperature gradient. In each slab 
compression is calculated by the Rankine-Hugoniot equations 
for the conservation of mass, momentum and energy through- 
out the shock front. In the case where shock and photoionization 
act on opposite sides the geometrical width of the nebula D is 
an input parameter, the diffuse radiation bridges the two sides. 
In this case, a few iterations are necessary to consistently obtain 
the physical conditions downstream. 

The ionizing radiation from an external source is charac- 
terized by its spectrum, which is calculated at 440 energies, 
from a few eV to KeV. Due to the radiative transfer, the radia- 
tion spectrum changes throughout the downstream slabs, each of 
them contributing to the optical depth. The calculations assume 
a steady state downstream. In addition to the radiation from the 
primary source, the effect of the diffuse radiation created by the 
gas emission (line and continuum) is also taken into account us- 
ing about 240 energies to calculate the spectrum. 

For each slab of gas, the ionic fractional abundances of each 
chemical element are obtained by solving the ionization equa- 
tions. These equations account for the ionization mechanisms 
(photoionization by the primary and diffuse radiation, and colli- 
sional ionization) and recombination mechanisms (radiative, di- 
electronic recombinations) as well as charge transfer effects. The 
ionization equations are coupled to the energy equation when 
collision processes dominate, and to the thermal balance if ra- 



diative processes dominate. This latter balances the heating of 
the gas due to the primary and diffuse radiations reaching the 
slab, and the cooling, due to recombinations and coUisional ex- 
citation of the ions followed by line emission, dust collisional 
ionization, and thermal bremsstrahlung. The coupled equations 
are solved for each slab, providing the physical conditions nec- 
essary for calculating the slab optical depth, as well as its line 
and continuum emissions. The slab contributions are integrated 
throughout the nebula. 

Dust grains are coupled to the gas across the shock front by 
the magnetic field, and are heated by radiation from the WD and 
collisionally by the gas to a maximum temperature which is a 
function of the shock velocity, of the chemical composition and 
radius, up to the evaporation temperature (Tdust ^ 1500 K). The 
grains throughout shock fronts and downstream might be de- 
structed by sputtering. Grain emission and absorption processes, 
which depend on wavelength and on grain characteristics, may 
therefore explain the complex structure of the light curves at dif- 
ferent wavelengths. 

To better understand the modelling of the SED, one should 
recall that the frequency of the bremsstrahlung peak depends on 
the shock velocity, while the radiation from an external source 
affects the SED between 10'"^ and 10'^ Hz. In a consistent way, 
the dust reradiation peak shifts to higher frequencies the higher 
the shock velocity, while the dust-to-gas ratio determines the ra- 
tio between the dust peak and the bremsstrahlung at a certain fre- 
quency. The reradiation of dust within the SS nebulae is, there- 
fore, observable in the IR only for relatively high dust-to-gas 
ratios, being the IR range dominated by radiation from the red 
giant (Contini 1997, fig. 4; Contini & Formiggini 2003, fig. 5). 

In the following, we adopt a more physical approach by mod- 
elling the SEDs at different epochs instead of light curves at dif- 
ferent wavelengths and translating magnitudes vs. wavelengths 
in fluxes (erg cm"^ s"') vs. frequencies. 



4. 1978-1986 active phase: the models 

Since the emission lines better constrain the physical properties 
of a nebula than the continuum, we start by modelling the com- 
posite line spectra. 



4. 1 . The line spectra 

Permitted and forbidden lines were observed along the whole 
outburst episode. Hack et al. (1988) reported that in quiescence 
only Ha and H/3 appeared as broad, relatively strong emissions, 
while the other Balmer lines were faint absorbed and severely 
blended with the neutral metallic lines. The FWOM ofH/3 varied 
from about 400 km s" ' at the beginning of the outburst to 600- 
700 km s^'from early 1978 to the middle 1981; then it remained 
at 600 kms^'until mid 1984, increasing to 750 kms^'in July 
1984; in November 1984 and 1985 it reached about 1200 km s"' . 
At the same epoch, high resolution UV spectra showed a strong 
and broad (4000 km s" ' ) Lya emission, never present at previous 
epochs. The origin of the broad Lya is discussed by Contini et 
al. (2007). 

Before 1985 the overall spectra did not show strong emission 
features except [SII] 4068+, [OI] 6300+ in the optical and OI 
1304, 1358, Silll] 1892, CIII] 1909, and Mgll 2800 in the UV 
(Mikolajewska et al. 1988). Emission lines of Hel 5876 were 
present during the whole outburst except in the beginning. The 
[OIII] 4363 and 5007 lines appeared only in November 1984, 
when the outburst was almost over. At the same epoch the UV 



Table 1. The models 
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no (cm ^) 


1.5e9 


2.e9 


5.e7 
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2.e7 
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70 
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2.e-3 
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3.e-3 


3.e-3 
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l.e-3 


l.e-3 


T.(K) 


1.5e5 


1.5e5 


1.5e5 


1.5e5 


3.5e4 


3.5e4 


3.5e4 
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0.5 


1 


6 
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0.01 


0.001 


0.001 
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4.0el4 


3.6el3 


l.el4 


l.el4 


5.el5 


3.el7 


3.el7 


dif 


4.e-4 


4.e-4 


4.e-4 


4.e-4 


4.e-7 


4.e-7 


4.e-7 



° radiation dominated (Ut^O); upper limit in radiation-bound models; ' by mass 



spectrum indicated that several high ionization emission lines 
(NV, CIV, SilV, Hell, OIII]) had appeared (Selvelli & Hack 
1985, Mikolajewska et al. 1988). The decrease in density at the 
end of the outburst was indicated by the ratio of FeII/[FeII] and 
CIII] 1906/1909. 

The UV line ratios suggested electron densities of 10^-10'" 
cm~^, characteristic of winds close to the stars. On the other 
hand, the forbidden optical lines indicated densities < 10^ cm"^ 
and were therefore emitted from the nebulae downstream of the 
expanding shock front. 

Since the structure and evolution of the nebulae are very 
complicated, we adopt three schematic types (Tables [I©: 
1) downstream of a shock with velocities of 600 - 1200 
km s 'corresponding to the reverse shock, 2) downstream of a 
shock expanding outwards with velocities < 150 km s 'which 
can be either radiation dominated, namely, the nebula is reached 
by the photoionizing flux from the WD, or shock dominated, i.e. 
the conditions in the nebula are those corresponding to a pure 
shock, and 3) downstream of a shock with velocities of ~ 40- 
60 km s^'and a high no (~ 10^ cm"^) propagating outwards the 
red giant star, which represents its own dust shells. As a matter 
of fact, the shock velocities of the shells are < 30 kms"', but 
the velocity increases when the shock propagates throughout the 
circumstellar environment with a negative density gradient. 

We can now explain the evolution of the spectra in the 1978- 
86 active phase. A complete grid of models was run and the mod- 
els which appear in Table[T]were selected from those which best 
fit the observed line and continuum spectra. 

In Fig.[3]we show the profile of the electron density, electron 
temperature and of the fractional abundances of the main ions, 
which explain the spectra. 

In particular, each diagram represents the profile of the phys- 
ical conditions in the gas downstream of a shock front. We have 
chosen three very different models. The diagrams correspond to 
shocked nebulae. The shock front is on the left. The diagram on 
the left of Fig. [3] and the middle one show the case in which ra- 
diation from the WD reaches the very shock front edge. In fact 
the diagram on the left represents the nebula downstream of a 
gas-Hdust shell ejected by the Mira. It is characterized by a high 
preshock density and a low Vj. The shell propagates towards the 
WD. 

The middle diagram represents the conditions downstream 
of the reverse shock, facing the WD. The shock velocity is 
roughly determined by the FWOM of the line profile, while the 
range of the preshock density is dictated by the line ratios. 

The diagram on the right of Fig. [3] represents the nebula 
downstream of the shock front propagating outwards the sys- 
tem. The X-axis scale is logarithmic and symmetric with respect 
to the middle of the nebula in order to show the conditions in 
both edges with similar high resolution, because, in this case, the 
shock and the radiation flux act on opposite edges. The physical 
conditions on the two sides of the nebula (the collisional domi- 



nated and the radiation dominated) are bridges by the secondary 
diffuse radiation from the hot slabs of gas. 

The middle diagram of Fig. [3] and the left side of the di- 
agram on the right show the characteristics of a shock domi- 
nated regime, namely, the temperature of the gas decreases from 
T - l.SlQr'iVsllQQkms^^Y immediately behind the shock- 
front. However, the recombination coefficients are low at high 
temperatures. Therefore, a long region of high temperature gas 
survives downstream, emitting, by the way. X-rays and coronal 
lines. The cooling rate is proportional to n^. When the temper- 
ature reaches about 10^ K the bulk of the UV and optical lines 
are emitted leading to a strong decrease of the temperature. The 
abrupt cooling is followed by the calculation code by reducing 
the geometrical thickness of the slabs (down to even 10^ cm) 
by (r,_i - Ti)/Ti-i < 10, where 7", and Tj^i represent the tem- 
peratures of two close slabs. After the drop, the temperature is 
maintained at about 10'* K by the primary and secondary radia- 
tions. Recall that photoionization alone hardly heats the gas to 
temperatures > 2-3 10"* K. In a pure shock case (U=0) the tem- 
peratures decreases rapidly below 1000 K and the neutral lines 
are thus strong. 

The situation is different in the shell (left diagram) with a 
relatively low temperature even close to the shock front. The 
radiative ionization rates prevail on the collisional ones which 
depend strongly on the temperature, so the picture is very similar 
to that of a pure photoionization regime. The density is very high 
and the optical thickness of the slabs is high enough to reduce 
the primary and secondary fluxes at a certain distance from the 
shock front. The gas therefore recombines. The case presented 
for the shell is matter-bound due to fragmentation at the shock 
front. 



4.1 .1 . Spectra from the reverse shock 

As mentioned before, in the epoch corresponding to the outburst 
the FWHM of the Balmer lines evolve from 400 to 1200 km s ' . 
These lines are emitted from the nebula downstream of the re- 
verse shock bounding the wind-wind collision region between 
the stars (Contini 1997, Angeloni et al. 2007a). Model m6oo is 
adopted to represent the high velocity shock: it shows relatively 
high Hfi , Hell and Hel line ratios (Table|2l). 

The FWHM of the broad line profiles suggests the range of 
the shock velocity Vj. High velocities are characteristic of the 
reverse shock between the stars, a high no is dictated by the range 
of densities in the star atmosphere. 

Recall that the temperature in the immediate post-shock re- 
gion is oc y2, so the temperatures are high (> 2.4 10^ K) for 
Vs> 400 kms"', recombination coefficients are low, and even 
compression (n/no) downstream is low close to the shock front, 
because the densities would correspond only to ~ 4 no (the adi- 
abatic jump). The radiation flux related to the H and He lines 
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Fig. 3. Left diagram: The profile of the electron temperature and electron density throughout the dusty shell for model mlsheii is 
shown in the top of the diagram. The fractional abundance of the most significant ions is shown in the bottom of the diagram. 
Middle: the same as in the left figure throughout the nebula downstream of the reverse shock corresponding tothe model (msoo). 
Symbols are the same as in the left diagram. Right: the same as in the diagrams on the left for the expanding nebula (ml^^p). The X 
axis units are cm. 



and the continuum can be as strong as to be observable only 
at temperatures < 10^ K. The densities play a dominant role in 
the cooling rate (oc n^, where n is the density) in order to re- 
cover those temperatures downstream of a strong shock. In fact, 
the density must be high enough (no > 10^ cm"^) to reduce the 
temperature to T< 10^ K within the nebula downstream of the 
reverse shock. This nebula therefore should have a geometrical 
thickness D smaller than the binary separation (10''' - 10'^ cm). 

This shows that Vj, no, and D are the main input parameters 
constraining the choice of the model. Moreover, a strong tranver- 
sal component of the magnetic field can reduce compression and 
a high dust-to-gas ratios may speed up the cooling rate of the gas 
through the collision heating of dust grains (Contini 2004). The 
ionization parameter and the star temperature aff'ect particularly 
the ionization structure of the gas. 

The radiation flux emitted from the gas downstream is cal- 
culated in each slab downstream and integrated throughout the 
nebula. The input parameters are even more constrained by com- 
paring the calculated flux Fc with the flux Fo observed at Earth 
throughout the continuum SED (Fig.|4|i. Recall that F^. r^=Fo d^, 
where r is the distance of the nebula from the system center and 
d is the distance to Earth. For instance, F^ = 2 10** erg cm"^ s^'at 
v=5 10'** Hz is calculated downstream of a shock front with 
Vs=600kms~'andno = 5 10^ cm"^. Forr~ 2 10'"* cm and adopt- 
ing d= 270 pc, Fo results ~ 10"^ erg cm"^ s"' , in agreement with 
the observations (Fig. IH. In this case the emission lines corre- 
sponding to the heavy elements are very weak because the low 
temperature (< 10^ K) zone in the nebula is narrow, while the 
Balmer and permitted He lines are high enough to be observed, 
due to a high ionization parameter U (> 1). 

The parameters are refined in order to obtain the best fit to 
observations in different epochs, considering that V, increases 
with time during the 1978-1986 active phase and therefore no 
must decrease. 



4.1 .2. Spectra from the shells 

At the earliest epochs (1978-1984)the models are constrained by 
the relative strength of OI 1304, 1358, Silll] 1892, CIII] 1909, 
and Mgll 2800. The spectra before November 1984 showed rel- 
atively strong low ionization lines which are explained by model 



ml,,,,,, (Tables [U and H with Vs=40 kms"'. The CIV/CIII] 
line ratio depends on the geometrical thickness D of the nebula 
namely, the larger D, the lower the CIV/CIII] ratio. Therefore, 
the observed CIV/CIII] < 1 suggests to represent the shells by 
matter-bound models relative to a quite extended nebula {D-A 
10'"* cm). 

The parameters describing the shells are chosen considering 
that dust shell characteristic velocities are ~ 30-50 km s"' and 
that the densities must be suitable to the giant atmosphere (> 
10^ cm"^). They are slightly modified in order to obtain line ra- 
tios fitting the observations. Notice that adopting low Vj (30-60 
kms"') and high densities, the temperatures downstream are < 
6 10"* K and the permitted UV lines emitted from the heavy ele- 
ments are high enough to be observed. 



4.1 .3. Spectra from the expanding shock 

The eclipse (Skopal et al. 1996) at the end of 1984 indicated 
that the system had rotated and the spectra were hence emitted 
from a different nebula, namely, the nebula downstream of the 
expanding shock. 

In November 1984, [OIII] and [NIII] in the optical domain, 
and several high ionization lines (e.g. CIV, SilV) appeared. The 
models best fitting both line and continuum spectra are shown in 
Table[3]where the line ratios are compared with the observations 
of Mikolajewska et al. (1988). 

The parameters corresponding to the expanding shock are 
dictated by the FWHM of the forbidden lines, by the densities 
that must avoid collisional deexcitation of the forbidden lines 
(e.g. [OIII] 5007), and by the consistent fit of the observed con- 
tinuum SED and of the Une ratios. 

In particular, the UV line ratios to Lya =10 calculated 
by model mlf^p are compared with the observations by 
Mikolajewska et al. (1988, Table© at JD 2446086 (20/1/1985). 
This model is shock dominated, as the spectrum is emitted from 
gas in the external side of the expanding nebula, facing the ob- 
server. Notice that a perfect fitting of the observations is not sen- 
sible due to the approximation of the models. Therefore we have 
focused on the line ratios which inequivocally describe a physi- 
cal situation. 
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Table 2. Calculated line ratios to H/3 



line 


^^ shell 


m2j,,f„ 


meoo 


mK„ 


NV 1240 


0.003 


0.23 


0.001 


4.15 


01 1304 


0.04 


0.02 


0.006 


0.006 


CII 1335+ 


0.8 


1.65 


0.15 


3. 


OI 1356 


0.42 


0.19 


0.06 


0.05 


OIV]+ 










SiIV1400 


0.34 


1.8 


0.02 


8. 


NIV] 1486 


0.03 


0.27 


0.002 


3. 


Sill 1531 


0.06 


0.03 


0.013 


0.02 


CIV 1550 


1.0 


7.2 


0.09 


14. 


Hell 1640 


3.8 


4.7 


1.9 


0.06 


OIII] 1663 


0.25 


0.9 


0.05 


2. 


NIII] 1744+ 


0.1 


0.4 


0.002 


1.3 


Sill 1814 


0.37 


0.2 


0.1 


0.08 


Silll] 1892 


0.7 


1. 


0.4 


2.5 


CIII] 1907+ 


1.2 


2. 


0.08 


5.6 


CII] 2326 


1.5 


1.2 


0.25 


3.2 


Mgll 2798 


1.1 


0.5 


3. 


0.5 


[Nelll] 3869+ 


0.008 


0.004 


0.001 


0.3 


[OIII] 4363 


0.008 


0.004 


- 


0.45 


HeI4471 


0.04 


0.02 


0.03 


0.63 


Hell 4686 


0.55 


0.65 


0.28 


0.005 


[OIII] 5007+ 


0.001 


0.002 


- 


1.47 


Hel 5876 


0.14 


0.08 


0.1 


1.12 


[OI] 6300+ 


0.002 


0.0004 


0.001 


0.25 


Bfl'' 


7.6e8 


1.9e8 


5.8e8 


0.75 



" shock dominated 
* inergcm"^ s"' 



Table 3. 1978-86: modelling the line ratios at JD 2446086 
(20/1/1985) 



line 


obS6086 


m2.,;,(,;; 


ml°,„ 


HI 1215 


10 


10 


10 


NV 1240 


1.6 


0.02 


0.8 


01 1304 


3.3 


0.02 


0.01 


CII 1335+ 


0.87 


0.13 


0.6 


OIV]+SiIV 1400 


0.83 


0.15 


1.3 


NIV] 1486 


0.19 


0.02 


0.43 


Sill 1531 


0.77 


0.002 


0.02 


CIV 1550 


1.8 


0.58 


2.7 


Hell 1640 


1.4 


0.4 


0.01 


OIII] 1663 


0.69 


0.007 


0.4 


NIII] 1744+ 


1.1 


0.03 


0.2 


Sill 1814 


1.23 


0.02 


0.02 


Silll] 1892 


0.89 


0.08 


0.5 


CIII] 1907+ 


1.2 


0.16 


1.1 



shock dominated 



The model is constrained by the CIV:CIII]:CII] line ratios. 
The shock dominated model mlg^-j, cannot explain the observed 
Hell 1640 line, which is better approximated by model ml shell- 
Also in the optical range model valsheii shows Hell 4686 and Hel 
5876 lines as intense as to be observed (Table|2]), while OI 1304 
and Sill 1814 are underestimated. In conclusion, also the shells 
contributed to the line spectra as well as to the continuum (see 
Sect. 4.2). 



the optical-UV range is zoomed in. The two bottom panels fo- 
cus, respectively, on the IR and optical-UV flux variations: in 
particular, the bottom left panel shows the time-evolution of the 
SED in the near-IR, characteristic of a significantly variable ab- 
sorption by dust. 

Model results are plotted over the data in the top panels of 
Fig. m The models are cross-checked by the line spectra pre- 
viously analysed (Sect.4.1) and constrained by the strong self- 
absorption in the radio range (top left panel). The high veloc- 
ity models mgoo and mi 200 are needed to explain the X-ray 
datum (top left panel). The models characterised by Vs> 600 
km s 'reproduce the Balmer line FWOM, but lead to UV-optical 
lines negligibly weak, so they were not included in Tables |2] and 

m 

Dust is collisionally heated by the gas to high temperatures in 
a strong shock velocity regime and sputtering reduces the grain 
sizes; the smaller the radius, the higher the temperature reached 
by the grains. Reprocessed radiation by dust downstream of the 
shock fronts appear in the top left diagram of Fig.|4] Model meoo 
can reproduce both the optical and radio continuum slopes. 

The dusty shells which fit the data in the optical-UV domain 
at maximum of activity (JD 2444937-2445720 = 28/11/1981- 
20/1/1984) contribute also to the CIV and CIII] lines in agree- 
ment with observations. Dust in the shells is mainly heated by 
radiation from the hot star: this was already found by previ- 
ous investigations (Munari et al. 1996, Angeloni et al. 2007b,c). 
Throughout expansion, the dust shells gradually fragmented by 
interaction with the environment, may reach the external edge 
of the binary system (> lO''* cm), becoming circumbinary and 
playing a primary role in the complex obscuration episodes of 
CH Cyg. 

At JD 2446212 (26/5/1985) and 2446776 (1 1/12/1986 - Fig. 
m right panels) at the end of the active phase a flat compo- 
nent appears in the near-UV SED. This is interpreted by model 
mlexp which represents the nebula downstream of the shock 
propagating outwards the system. Model val^xp is shock dom- 
inated (U=0) and explains most of the line ratios observed 
by Mikolajewska et al. (1988) at JD 2446086 (20/1/1985 - 
Table O. The slopes at JD 2446212 (26/5/1985) and 2446776 
(1 1/12/1986) can also be fitted by the same ml„p model in the 
radiation dominated case, adopting a relatively low U. However, 
in this case, the calculated line ratios are lower than observed. 
Recall that the line spectrum and the continuum were observed 
at different, although close, days. This is a further evidence that 
if the medium surrounding the system is clumpy enough, the 
radiation flux throughout its path from the WD to the external 
shock can be flickering. 

During the brightness minimum in 1985, CH Cyg showed 
X-ray emission (Leahy & Taylor 1987). Simultaneous radio ob- 
servations with the Very Large Array (VLA) revealed a radio 
outburst and expanding jets (Taylor et al. 1986). The X-ray da- 
tum in Fig.|4]was taken at the end of 1984, therefore constraining 
the models at that time. For earlier times, when no data are avail- 
able, Leahy & Taylor suggested a flux lower by a factor of ~80. 
We agree with them that the very large column density of the 
downstream nebulae led to large absorption in the 0.05 -1.5 ke V 
range. 



4.2. The continuum SED 

In Fig. |4] we present the modelling of the observations (Table 
lA.lb at some significant days during the 1978-86 active phase. 
The top left panel of Fig. |4] displays the data on a large fre- 
quency range (from radio to X-ray), while in the top right one distance to Earth. 



4.3. Radii of shells and nebulae 

Comparing the bremsstrahlung calculated at the nebula with the 
data observed at Earth we obtain the factors 77 = r^/d-, where r 
is the distance of the nebula from the system centre and d is the 
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The reverse shock (model mgoo) was hence located at a ra- 
dius of 2 lO'"^ cm. The top right panel of Fig. |4]points out that 
the bremsstrahlung from the shell in the optical-UV range, rep- 
resented by model m2si„ii, can explain the data at two differ- 
ent days, JD 2443877 (3/1/1979) and JD 2444937 (28/1 1/1981). 
The calculated radii are 2.88 lO'-' and 9 10'"* cm, respectively, to 
which would correspond a shell expansion of ~ 6 10'^ cm in 2.9 
years. This, in turn, would correspond to a velocity < 7 kms ', 
much lower than the shock velocity (Vs=60 km s"'), unless the 
two shells were ejected at different times. The latter hypothe- 
sis seems more realistic, as the pulsation frequency of the cool 
component is actually ~756d. 

The expanding shock reached a radius of 8 10'^ cm. The 
geometrical thickness A of the cool region (T< 10"* K) within 
the downstream nebula is constrained by the relatively strong 
self absorption of free-free radiation in the radio. The opti- 
cal thickness r is oc n^ A, where the density n (< 10^ cm""' 
considering compression downstream) is constrained by the 
slope of the bremsstrahlung: a larger n would lead to a steeper 
bremsstrahlung in the UV-optical range. This allows us to give a 
lower limit to the geometrical thickness of the cool region ; A > 
2 10'^ cm. 
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Fig. 4. Modelling the active phase 1978-86. Top pannels: com- 
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Fig. 5. Modelling of the quiescent phase 1987-90. The figure is 
structured as Fig.|4l Top panels: Solid line: bb (2000 K); short- 
dashed lines: m3exp', long-dashed lines: m2„p; dot-dashed lines 
: model val^^p. The models refer to Table [1] 



5. 1987-1990 quiescent phase: the obscuration 
effect 

5. 1 . The line spectra 

Taranova & Yudin (1992) stated that the optical spectra taken 
in 1988 did not reveal any emission lines (Mikolajewski et al. 
1990) or revealed only very week Ha and H/3 lines (Andrillat 
1988, Bopp 1990). In 1989, before the next outburst, an emis- 
sion line spectrum with Fell, [Fell], [SII] 4068, and faint one- 
component Hy6 was lying upon numerous absorption lines. The 
spectrum can be roughly explained by the expanding nebula rep- 
resented by model ml„p with Vs=150 kms"'. At JD 2448122 
(18/8/1990), Tomov & Mikolajewski (1992) observed numerous 
emission lines from the shell which can be explained by model 

ni2.5/,e;;. 



5.2. The continuum SED 

In Fig. |5] we present the SEDs observed during the quiescent 
period 1987-1990 between the two active phases 1978-86 and 
1991-95. Fig. |5]is structured as Fig. |4] In the top right panel 
the V data are fitted by the bb flux from the cool star, while 
the B and U data showed a flat slope that can be reproduced by 
models ml(,^p and m2j,vp, corresponding to dififerent preshock 
densities (Table[Tli. So we constrain the models by the line ratios, 
recalling (Sect. 5.1) that the optical spectra observed in 1988 did 
not reveal any emission lines or very week Ha and H/3 . This 
excludes the shock dominated model mlexp that was used to fit 
the flat slopes at the end of the 1978-86 active phase and the line 
spectrum at JD 2446086 (20/1/1985). 

The models which best fit the slopes at JD 2448136 
(1/9/1990) and JD 2447339 (26/6/1988), m2„p and m3„p, re- 
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spectively, are both shock dominated with Vs< 100 km s 'and 
relatively high no, which lead to negligible line emission down- 
stream. They indicate that the outflowing shock is propagating 
through a denser medium, perhaps a residue of matter ejected 
during previous eruptions. 

Finally, let us compare the near-IR fluxes during the 1978-86 
phase with the ones during the 1 987-90 phase (bottom left panels 
of Figs.|4]and|5] respectively): while in the former case the slopes 
showed a dust absorption clearly wavelength-dependent, dur- 
ing the quiescent phase minimum the IR flux decreased equally 
throughout the whole near-IR range, suggesting an obscuration 
effect by an external shell. 

6. 1991-1995 active phase: the black body radiation 
from the WD 

A new active phase began in December 1991 with a large flux 
rise in U, a smaller one in B and a negligible one in V. This trend 
was sustained until March 1992. Then, from the beginning of 
1993 the U band rose again until summer 1994. The eventual 
drop in summer 1995 informed that the outburst was running 
out. 

It is worth noting that during this 1991-95 high state, the bb 
continuum flux from the hot star became gradually detectable 
(right panel of Fig.|6l). As a matter of fact, the appearance of a hot 
continuum source with T(.//= 150,000 K was taken as the indica- 
tion of the beginning of the outburst by Skopal et al. (1996). To 
obtain an initial hint about the system during this active phase, 
we compare the evolution of the slopes in the U-B-V during the 
1978-87 and 1991-95 active phases in Fig.|6] An impressive dif- 
ference can be noticed comparing the right and left panels: in 
particular, the left one shows that in the first active phase the 
trend in the UB range cannot be explained by the bb flux from 
the hot star which vice-versa is evident (right diagram) during 
the 1991-95 outburst. The direct contribution of the hot star to 
the U band is shown at the 1991-95 active phase maximum (JD 
2449598 = 2/9/1994). In Fig.|6l the low frequency tails of the 
black body curves correspond to T»= 100,000 K (dashed line) 
and T» =35,000 K (dotted line). These temperatures represent the 
maximum and minimum one, respectively, during and after the 
outburst (Corradi et al. 2001). The slopes of both bb fluxes in the 
UB frequency domain are slightly different than observed, and 
are even less fitting the data before (thin lines) and after maxi- 
mum (thick lines). We suggest that the WD was hidden from a 
direct view in the 1978-86 active phase, while it became visible 
during the 1991-95 maximum. 

The two deep minima in U recorded at the end of September 
1992 and in October 1994 lasted 81 and 127 days, respectively. 
They were attributed to the eclipse of the active component by 
the cool one in the light of the triple star model (Skopal et al 
1996), but they are hardly explained by the binary interpreta- 
tion which ascribes the 756d periodicity to a non-radial pulsa- 
tion of the M giant. Fig. |6] (bottom panel) shows that the eclipses 
of September 1992 just before JD 2448900 (4/10/1992) and of 
October 1994 involved only the B and U bands, the WD and the 
shocked nebula being eclipsed most probably by dust shells. The 
V band in fact was not involved as it recovered its range of mag- 
nitudes that were shifted during the bursts by the contribution of 
the nebulae. Similarly, the relatively weak maxima immediately 
after the "eclipses" were due to emission from the shocked neb- 
ulae. We suggest that the minima were due to absorption by the 
dusty shells that are ejected in tune with the pulsation of the red 
giant. Recall that the dusty shells are highly fragmented there- 
fore this effect cannot be observed periodically. 



6. 1 . The line spectrum 

At JD 2448094 (21/7/1990), a dramatic change led to broad hy- 
drogen absorption profiles and at JD 2448095 (22/7/1990) CH 
Cyg displayed faint emission in Ha with corresponding veloci- 
ties around -67 kms"'. Then, in August 1991, the appearance of 
absorption terminal radial velocities of - 1 000 and -900 km s ' for 
Ha and H/3 , respectively, might have been a consequence of the 
rotation of the system. This would have allowed us to see the 
reverse shock in front of the WD. Eventually, the absorptions 
turned into emissions with velocities up to 1800 for Ha and 1600 
kms^'for H/3 : these velocities disappeared at the beginning of 
flie 1992. From JD 2448883 (17/9/1992) the H/3 emission in- 
creased, whereas no absorption lines were observed. 
High velocities had been observed also during the previous ac- 
tive phase (Sect. 4) and there explained by the reverse shock 
model. 

The spectral region around [OIII] 5007 showed complicated 
structures consistent with more components. During the outburst 
(1992 onwards) strong emission lines of Fell, [Fell], [SII], and H 
were observed, while [NeIII]3869+ and [Felll] 4659-4670 were 
still faint, increasing later. Also, [OIII]4363 and [OIII] 5007+ 
were relatively faint and varied rapidly during the active phase. 
The overall scenario, with multiple emission line components 
and structured line profiles, may be qualitatively understood by 
Richtmyer-Meschkov (R-M) and Kelvin-Helmoltz (K-H) insta- 
bilities (Graham & Zhang 2000). R-M is an interfacial insta- 
bility between two fluids of different density driven by shock 
waves. The K-H instability leads to the consequent fragmenta- 
tion of matter at the shock front (Contini & Formiggini 2001). 
The dumpiness of the nebula might then contribute to flickering 
episodes. 

On September 1992, the M-giant was even more apparent 
due to the eclipse. The red part of the spectrum was dominated 
by the M-giant molecular spectrum, overlayed by strong Ha , 
[OI]6300, and Hel 5876 emissions. There were simultaneous 
emission lines of [SII], [Nelll], and [Felll]; notice that the Hell 
4686 and 1640 lines were not present during the outburst. We 
suggest that the line spectrum is most probably emitted from 
a disrupted shock dominated expanding nebula (ml^vp) charac- 
terised by low Hell lines (Table|2]). 

In May 1993 the UV emission line spectrum showed Lya , 
OIV], SilV 1400, CIV 1550, OIII] 1666, Silll] 1892, CIII] 1909, 
CII] 2325, [Oil] 2470. These lines are reproduced by the model 
for the expanding nebulae ml,.^^. Moreover, numerous Fell mul- 
tiplets appeared, together with Mgll 2796, 2803 doublet, and 
sharp emission lines, Arlll 2187, 2191. In August 1994, a sharp 
absorption cutting the emission line profiles at Vs=150 kms"' 
most likely shows that the expanding nebula is facing the ob- 
server In fact, the strongest emission came from the line core 
while a very broad and smooth absorption up to - 1 600 km s ' was 
seen. 



6.2. The continuum 

The SEDs of the continuum are shown in Fig. |7] diagrams. The 
models which appear in the left diagrams, covering the whole 
spectral range, were selected by the consistent fit of the data in 
the optical-UV range (zoom in the middle diagram), of the radio 
and of the X-ray data. 

We start by presenting the modelling in the optical-UV 
range because it is more constraining. Around maximum at JD 
2449598 (2/9/1994), the flux SED corresponding to the sum of 
the expanding shock (ml^^^^) bremsstrahlung and the bb flux 
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et al. (2000) claim that flickering cannot be observed during the 
eclipse. 

Similarly, the flux SED which best reproduces the trend of 
the data at JD 2449227 (27/8/1993) results from the sum of 
model mi2()() with the WD bb flux: recall that broad lines were 
still emitted in 1993. In May 1995 we observed only the expand- 
ing shock corresponding to model mlg^p. 

Fig El middle diagram shows that the bb flux from the WD 
grew by a factor of 2.75 between JD 2449227 (27/8/1993) and 
JD 2449598 (2/9/1994), corresponding to an increase of the WD 
temperature from ~ 1 15,000 K to 150,000 K. 

We can now examine the data on a large scale, from radio to 
X-ray. The radio slope between 1992 and 1993 as observed by 
MERLIN, VLA, and JCMT (Skopal et al 1996) is once again 
explained by bremsstrahlung with a strong self-absorption at 
longer wavelengths. The best fit is obtained with a model cor- 
responding to Vs=1200 kms 'in the reverse shock. This model 
also explains the data in the optical-UV range and the hard X- 
ray. Models mlexp and mlexp give a negligible contribution in 
the radio. Skopal et al. (1996) proposed a velocity of ~ 310 
kms" 'in order to match the contemporary angular radius ob- 
served by MERLIN: this discrepancy is easily explained as we 
deal with disrupted nebulae downstream of head-on shocks in a 
highly turbulent regime. 

The soft X-ray datum by ASCA was modelled by Ezuka et 
al. (1998) as two optically thin thermal plasma emissions. We 
find (Fig. |7] left diagram) that the soft X-ray is emitted down- 
stream of the expanding shock, while the hard X-ray would cor- 
respond to the reverse shock between the stars with a velocity of 
1200 kms'. 

An impressive systematic evolution of the JKLM band flux 
appears in the right diagram of Fig. [T] indicating that dust is 
involved in the obscuration process. 



from the hot star (T»= 150,000 K) reproduces the slope of the 
U-B-V SED. At JD 2449583 (18/8/1994) the bb component was 
not visible. This may be explained by the fragmented nature of 
the shells which screen the flux from the WD, confirming Eyres 



Finally, by exploiting the 77 factors defined in Sect. 4.3 we 
can calculate for the reverse shock a radius of 4 10'^ cm and for 
the expanding nebulae which best fit the data a corresponding 
radius of 10'^-3 10'^ cm. Similar dimensions of the system were 
found by HST imaging in the 1998-2001 active phase by Eyres 
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et al (2002) after about 5 years. Adopting en expanding velocity 
< 100 kms~'the radius increment is relatively small (< 2 10'^ 
cm). 



7. 1996-1997 transition phase: the contribution of 
radio synchrotron radiation 

During the quiescent phase 1996-97 some interesting radio ob- 
servations were taken between March and October 1997 by 
Karovska et al. (1998). The data (Fig. |8] top left), along with 
an intriguing flux variability, followed the trend of an absorbed 
bremsstrahlung as was found in previous phases adopting a high 
velocity model (Vs=1200 kms"'). A particular investigation of 
the radio range is shown in Fig. [8] bottom left diagram. The 
slope of the data at 8 and 43 GHz observed between May and 
August 1997 could be fitted by synchrotron emission created at 
the shock front by the Fermi mechanism in the non-relativistic 
case (Bell 1978) with a characteristic spectral index a - -0.35. 
This is not surprising because non-thermal emission was also in- 
voked by Crocker et al. (2001) in order to explain the radio data 
between 1985 and 1999, adopting a magnetic field of the order 
of 1-3 mG. Interestingly, we have reached the same values of the 
magnetic field (Bq = 1-3 mG, Table[T]i by consistent modelling of 
line and continuum spectra in the different phases. Nonetheless, 
the datum at 22 GHz is highly underpredicted by the synchrotron 
radiation and may refer to a completely different emission pro- 
cess, e.g. bremsstrahlung from the downstream nebula. A high 
velocity model is in fact requested by the optical-UV data, as 
presented in Fig. [8] (top diagrams). 

Regarding the infrared range, it is worth reminding some re- 
sults of 1997 ISO observations, unfortunately still no fully ex- 
ploited. Schild et al. (1999), who did a preliminary analysis of 
the spectra, found no strong emission lines but [OI]63 /im, and 
strong OH and weak H2O molecular emission on top of a rapidly 
decreasing continuum. This continuum can be fitted by a single 
dust shell model (T^ ~ lOOOK, Fig. [8]), in contrast with other 
SSs which seems to display multiple dust shells with character- 
istic radius and temperatures (Angeloni et al. 2007c). 
The dust chemistry was quite surprising, since weak PAHs (at 
6.3 jum) coexisted with dominant amorphous silicate (at 9.7 and 
18 //m). Unfortunately, a comparative study of dust spectroscopy 
capable of investigating the grain formation and destruction pro- 
cesses throughout both the active and quiescent phases is still 
not possible for lacking of data (see Angeloni et al. 2007a for 
a comparison of emission lines and silicate dust features, from 
ISO spectra, in a various sample of SSs including CH Cyg). ISO- 
SWS observed the system twice during this quiescent phase, and 
in Fig.|9]we show a comparison of the two recorded spectra. As 
one can notice, the maximum flux variation seems to be in corre- 
spondence of the strong silicate emissions, while the underlying 
dust continuum appears not so sensitive to the environment vari- 
ation due to the giant pulsations. How this can be related to the 
actual dust formation and processing is far from understood. 

8. 1998-2001 active phase: jets and X-rays 

8.1. The line spectra 

As noticed by Eyres et al. (2002) in the HST imaging, the very 
inner regions of CH Cyg contributed mainly to continuum, and 
very little to the line emission. In our interpretation this is not 
surprising since the nebulae emitting the line spectra would cor- 
respond to the expanding nebulae that are moving towards the 
outside of the system. 
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Fig. 9. The two ISO-SWS spectra. Top panel: the spectra dis- 
played in arbitrary units. It is noticeable the variations in the 9.7 
/vm silicate profile. Bottom panel: the relative ratio of the fluxes, 
showing that the main variations seems to be due to the silicate 
emissions, and not to the underlying dust continuum. 
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Balmer lines Ha , H/3 , Hy, all double peaked, lay upon a 
broad feature 1200 km s"' wide confirming that a fast reverse 
shock of Vs=1200 km s"' was present within a small radius < 
10'^ cm. The fast shock can explain also the X-ray Chandra ob- 
servations of 2001 at ~ 0.26 and 2 keV reported by Galloway 
& Sokolowski (2004). Indeed, the temperature in the immedi- 
ate post shock region corresponding to that Vj, is 2.16 10^ K, 
decreasing downstream to ~ 1.5 10^. The cooling rate is not so 
strong because recombination coefficients are low at such high 
temperatures. 

The HST images by Eyres et al. (2002) revealed strong 
[OIII], which can be explained by the expanding nebulae miexp 
modelled by fitting the SEDs. 

In September 1998 the OI, Hel, Hell, Fell, and Till were 
probably emitted from the shells. [OIII] 5007 was also evident, 
while [OIII] 4363 absent both in and out the eclipse. We have 
run some further models corresponding to expanding nebulae 
(MEfip and MEl^p) which, in agreement with the observations, 
predict relatively low [OIII] 4363 (Table |4]l. The input parame- 
ters are shown in the bottom of Table |4] The models are shock 
dominated (U=0). The dust-to-gas ratios adopted in these mod- 
els are 4 10""* by mass. 

Corradi et al. (2001) indicated strong [Oil] and [Nil] Unes 
in the 1999 images. These strong low ionization lines are char- 
acteristic of shock dominated models, e.g. MEpvp and MElpvp, 
which might occur from the collision of the expanding jets with 
the surrounding ISM. The radius of these shocked nebulae are 
calculated by fitting the continuum SED. The rj factors lead to a 
radius > 3 10'* cm. 



8.2. The continuum SED 

In Fig. [To] we present the modelling of the 1998-2001 active 
phase. Some significant UBV data are taken from Eyres et al 
(2002, Fig. 2). Fig. [To] (left diagram) shows that the best fit of 
the UBV slopes are obtained by model mle^p in December 1998 
during the maximum, while at later epochs the slopes are repro- 
duced by model mi2oo- This can be explained considering that 
at maximum the relative position of the different system com- 
ponents (stars, nebulae, shells) is such that we cannot see the 
reverse shock, which appears only after some months. 

The fitting of the models to the data in Fig.[TO](right diagram) 
leads to similar factors 77 for the two models. However, we know 
that the radius of the reverse shock is much smaller than that of 
the expanding shock by at least a factor of 100. This question is 
resolved invoking the filling factor as a key parameter. The radii 
of the expanding shell at maximum of the active phase and of 
the reverse shock at later epochs result of 1.3 IQ^'ljf cm and of 
1.4 lO^^ljf cm, respectively, leading to corresponding filling fac- 
tors of ~ 0.001 and 0.1, to recover more sensible radii (cfr. Sect. 
4.3). Small jf are suitable to the fragmented character of matter 
throughout the expanding shock deriving from dynamical insta- 
bilities e.g. R-M and K-H (Sect. 6.1 and Contini & Formiggini 
2001). 

Fig. [To] (left diagram) shows also that the djg ratio adopted 
by model ME^ ^^ is just a lower limit, because the flux cor- 
responding to reradiation by dust is constrained by the SED, 
definitively dominated by bb from the cool star in that wave- 
lenght range. The low velocities of these models indicate that the 
expanding nebulae swept up a large amount of interstellar mat- 
ter, whose dig actually corresponds to djg > 4 lO""'. This matter 
probably existed before the start of this outburst, and might be 
considered as a sort of remnants from the former bursts which 
have then propagated throughout the ISM. 



Active phase 1998-2000 (JD 2450000+) 




Fig. 10. The 1998-2001 active phase. Left diagram: modelling 
the SED from radio to X-ray. Right diagram: the optical fre- 
quency range. Long-short-dashed lines: mnoo; dash-dotted lines: 
va2exp\ dotted lines: MEexp', solid line: bb (2000 K). Large black 
triangle: data from Karovska et al (2007). The vertical thin lines 
show the range of strong absorption in the X-ray domain. 

Table 4. Line ratios to H/? in the extended nebulae (jets) 



line 


ME„p 


ME1„„ 


[Oil] 3727+ 


40. 


19. 


[OIII] 4363 


2.6 


2.18 


Hell 4686 


0.016 


0.027 


[OIII] 5007+ 


35.3 


25.7 


Hel 5876 


04 


0.58 


[OI] 6300+ 


0.37 


0.40 


[Nil] 6584+ 


8.9 


6.42 


H0(ergcm-2s-') 


0.0013 


0.0025 


Vs(kms-') 


70 


100 


no (cm"'') 


5.e3 


5.e3 


Bo (gauss) 


l.e-3 


l.e-3 


D(cm) 


l.el5 


l.el5 



9. From 2001 up to date 

Since 2001, CH Cyg persists in a quiescent phase recovering 
from the latest high state, at a rather bright magnitude (U~10 
or less, V=7-8, B=8.7-9.4). During the last months it showed 
a ~ 2mag decline in the optical bands (2006 June-December, 
Skopal et al. 2007) and a most significant, anomalous, IR fading 
(Taranova & Shenavrin 2007). 

In April 2004, high-resolution optical data by Yoo (2006) 
detected double-peaked emission components for the Via and 
Wp profiles, with an interesting V/R <1, where V and R stand 
for the intensities of the blue and red wings, respectively. This 
has been taken as a proof, along with the vanishing of the 
blue continuum, that at that time an accretion disk was akeady 
present, although highly processed. Moreover, the multiple- 
peaked profiles of the emission lines led to suggest a strong non- 
homogeneous envelope, result of the collision between some 
outflowing matter and the circumstellar shells. 

A few months later, in October, new observation revealed 
(Yoo 2007) a rapidly developing framework, with single-peaked 
Hot and HyS emissions and no absorption component at their line 
centres. The V/R ratio became << 1, indicating that the accre- 
tion disk might have been destroyed within the April-October 
time interval. Variability in the Ho- blue component intensity 
has been explained as the actual expansion effect of some jets. 
Moreover, the blue continuum became so dim as to unveil the 
M giant even at short wavelengths, while a contemporary en- 
hancement of the emission lines allowed to detect conspicuous 
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Quiescent phase 2000-,.. (JD2450000 + ) 
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Fig. 11. Red line : the SED at maximum of U; green line : 
the SED at maximum of V; blue line ; the SED towards min- 
imum. Short dashed line: model mlexp', short-dash-dotted line: 

wl mlfvp + w2 bb2()oo, where wl = 10"'^'^ and w2=10"'^"* are 
the relative weights, respectively; long-dash-dotted: wl ml„p 
+ w2 bb2ooo where wl^lO"'-'-'^ and w2=10"'^° are the rela- 
tive weights, respectively; dotted: wl mlexp + w2 m3exp where 
wl = 10 '^^ and w2=10"'**^ are the relative weights, respec- 
tively. 



[OIII] 4959 and [OIII] 5007, along with broad emission compo- 
nents of [Nelll] lines. This emission might have originated in an 
extended, expanding shell constructed by the outflowing matter 
from the WD. 

In the X-ray, Mukay et al. (2006), analysing Suzaku obser- 
vations of 2006 January and May, found CH Cyg in a soft X-ray 
bright, hard X-ray faint state. While the soft X-ray flux level was 
similar to the 1994 ASCA observation, the hard one was much 
lower, confirming that the hard component is more variable than 
the soft component (Ezuka et al. 1997), which in turn may be 
spatially distinct. Assuming a binary period of ~14.5 years, the 
entire binary should be extended <0.4" (lOOAU at ~270 pc), ex- 
plaining why Chandra was not able to resolve the system and, 
at the same time, in agreement with an extended X-ray emis- 
sion region, comparable with the size of the binary, as expected 
in the colliding wind interpretation demonstrated by consistent 
modelling. 

The hard X-ray component, more directly sensitive to the 
variable symbiotic activity, is the manifestation of those non- 
thermal emission processes (such as synchrotron) which once 
again prove the key-role played by shocks in accelerating rela- 
tivistic particles responsible for the bulk of emission at hv >10 
(De Becker et al. 2007). 

In Fig.[TT|we present the modelling of the UBV and JHKLM 
data presented by Skopal et al. (2007) and Taranova & Shenavrin 
(2007), respectively. The models which best fit the slopes of the 
data are the combination of those used in previous epochs (see 
caption of Fig. fTTT i. Three main slopes can be noticed. That at 



maximum of U (red) shows a contribution of two nebulae with 
different densities. The curve with the highest V (green) have 
a conspicuous contribution from the giant bb radiation. This 
contribution decreases and even disappears at minimum (blue 
curve), indicating that the whole binary system is obscured by 
the expanding nebula, as was found in previous epochs. 



10. Discussion and concluding remarks 

The complex nature of the SS CH Cygni has been investigated 
by modelling the line and continuum spectra at different epochs, 
from the impressive outburst of 1978 up to date. Our results give 
an approximated picture of the shocked nebulae throughout the 
system because our approach implies a plane parallel geometry 
for the different shock fronts. This is constrained by the compos- 
ite models which account for both photoionization from the WD 
and shocks created by the collision of the winds. However, the 
plane-parallel geometry can be justified for relative large radii 
and in particular situations characteristic of SS. 

The present analysis has revealed that similar developing 
trends may be recognised throughout the active phases and 
throughout the quiescent phases. In the light of a colliding- wind 
scenario, we have specifically focused on the role of the nebu- 
lae created by collision of the winds from the hot and the cool 
stars, in order to build up a comprehensive view of the system 
as a whole, namely to self-consistently interpret the variable line 
spectra, the radio emission, the IR photometric properties, the 
fluxes in the UBV bands and the X-ray observations. The most 
important result of our method that can be applied to all the SSs, 
consists in recognizing the different sources of radiation (stars, 
nebulae, dust) throughout the continuum SED. 

We have found that bremsstrahlung from the nebula down- 
stream of the shock front facing the WD between the stars ex- 
plains the continuum slopes observed during the active states. 
These fast shocks seem to be always present during the out- 
bursts, proving that the wind collision is playing a primary role 
in the spectral properties of the symbiotic phenomenon (e.g. ra- 
dio and X-ray wavelength region). The circumstellar dusty shells 
define the optic-UV SED , also contributing to the line spectra 
during the active phases. 

The WD may be related to flickering during the active 
phases, either by appearing directly (1991-95) or by photoion- 
ising the shells throughout a disrupted medium (1978-87). 
Towards the end of the active phases, the expanding shock prop- 
agating outward the system is visible both in the continuum SED 
and in the line spectra. 

During the quiescent phases the circumbinary expanding 
shocks dominate. We have found that these shocks propagate 
throughout a relatively dense medium which may be interpreted 
as matter ejected in previous outbursts. This leads to a rela- 
tive strong synchrotron radiation from the shock front at some 
epochs. 

The jets observed by Corradi et al (2001) imply the presence 
of an accretion disk. The modelling of the low-ionisation lines 
(e.g. [Nil, [Oil]), that are strong in their images, confirms that 
the jets expanding in the ISM, have reached a radius > 0.01 pc 
from the SS. 

We have shown, by consistent modelling, that soft and hard 
X-rays are emitted downstream of shock fronts corresponding 
to quite different velocities, namely, the expanding shock (Vs~ 
100-150 km s"') propagating outwards the binary system would 
explain the soft X-ray emission, while the hard X-rays corre- 
spond to the high velocity shock between the stars, in agreement 
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with the recent results discussed by Karovska et al. (2007) on 
the basis of radio, optical and X-ray imaging. 

Dust grains within the nebulae are collisionally heated by 
the gas in a strong shock velocity regime (> 600 km s '), while, 
in the dusty shells, grains are mainly heated by radiation from 
the hot star The dust expanding disrupted shells might reach, in 
some cases, the external edge of the symbiotic system (r> lO''* 
cm), becoming circumbinary and therefore contributing to the 
complex obscuration episodes of CH Cyg. 

It should be emphasized that the interpretation of CH Cyg 
spectral evolution is further complicated by the relatively short 
orbital period (~14.5 years): in fact, the orbital motion combined 
with the unpredictable recurrent activity leads to different system 
configuration as seen from Earth. We have tried to explain them 
by combining the radiation fluxes from the several emitting con- 
tributors with different, variable weights. 

For the sake of clarity, we schematically summarise the re- 
sults of our modelling: 



1978-1987 active phase 



kms 



1200 



- Emission lines with FWHM of 400 
km s 'correspond to the reverse shock; 

- before 1984 November, low-ionisation lines come from the 
shells; 

- after 1994, higher level lines are emitted both by the expand- 
ing nebula and the dusty shells; 

- at days JD2446212 and JD2446776, UV and optical lines 
from the expanding nebulae appear; 

- at the end of the active phase, shock dominated models ex- 
plain both line and continuum spectra. 

1 988- 1 990 quiescent piiase 

- The expanding nebula propagates throughout matter ejected 
in previous eruptions; 

- in the IR the system is obscured independently from wave- 
length. 



1991-1995 active piiase 

- In 1991, at the onset of the new active state, the lines dis- 
play FWHM< 100 km s ', indicating a disrupted expanding 
nebula; 

- in August 1991, broad Ha and Hfi are seen in absorption, 
suggesting that the reverse shock is absorbing the WD flux; 

- by 1992 absorption lines disappear and M-giant molecular 
bands are observed. In May 1993, UV lines are emitted from 
an expanding nebula with a radius of 10'* - 10'^ cm. 



1996-1997 quiescent piiase 

- The radio flux shows a synchrotron component created by 
the Fermi mechanism at the shock front. 



1998-2001 active ptiase 

- . A fast reverse shock of 1200 km s 'may explain the X-ray 
Chandra observations in 2001; 

- the HST images by Eyres et al. (2002) reveal strong [OIII], 
which can be explained by the expanding nebulae in agree- 
ment with the fit of the SED; 



- in 1998 Sept the OI, Hel, Hell, Fell, and Till fines are emit- 
ted from the shells. Moreover, [OIII] 5007 is evident and 
[OIII] 4363 is absent both in and out of the eclipse. These 
lines come from shock-dominated expanding nebulae with 
relatively low velocities (~70-100 km s"') and relatively low 
densities (no~ 5000 cm"^). In fact the corresponding [OIII] 
4363 lines are relatively low; 

- the 1999 images of Corradi et al. (2001) indicate strong [Oil] 
and [Nil] lines: they are characteristic of a shock dominated 
nebula with a relative low no=5 lO"* cm"-'. 

2001 -up to date 

- three different slopes are observed in the optical-IR range. 
They are modelled by the combination of different expand- 
ing models with the bb flux from the cool component star 

- the decline observed by Taranova & Shenavrin in the IR is 
consistent with the decline in the optical-UV range and is 
explained by the eclipse of the whole system. 
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Table A.l. Observation references for the 1978-86 active phase. Table A.2. Observation references for the 1987-90 quiescent 

phase. 
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Observation 


References 




date 




radio 


1989 Mar 7-12 


AltenhofF etal. (1994) 




1988 Oct 14 


Ivisonetal. (1991) 


IR 


1988-90 


Taranova & Yudin (1992) 


optical 


1988-90 


Taranova & Yudin (1992) 


UV 


- 


- 


X 


- 


- 



Table A.3. Observation references for the 1991-95 active phase. 



Observation 
date 



References 



radio 1 992 Jul 24 Skopal et al. ( 1 996) 
1993 Jun 17 - Aug 19 Skopal et al. (1996) 

1993 Skopal etal (1996) 

IR 1991-95 Munari et al. (1996) 

optical 1989-95 Skopal et al. (1996) 

1989-94 Skopal etal (1996) 

1991-95 Munari et al. (1996) 

UV 1 989-95 Skopal et al. ( 1 996) 

1991-95 Munari etal. (1996) 

1 992 Nov 24 Skopal et al. ( 1 996) 

1 993 May 27 Skopal et al. ( 1 996) 
X-rays 1994 Oct 19 Ezukaet al. (1997) 



Table A.4. 

phase. 



Observation references for the 1996-97 quiescent 



Observation 
date 



References 



radio 


1997 Mar 8 -Aug 31 


Karovska etal. (1998) 


IR 


1997 May 10 


Schild etal. (1999) 


optical 


1995 Dec 18- 1997 Mar 11 


Skopal etal. (1997) 




1996-97 


Corradi et al. (2001) 


UV 


1997 Jul 10 


Corradi et al. (2001) 


X 


- 


- 



Table A.5. Observation references for the 1998-2000 active 
phase. 





Observation 


References 




date 




radio 


1999 Sep 26 


Eyres et al. (2002) 


IR 


1999-2003 


Taranova & Shenavrin (2004) 


optical 


1998 Sep 5 -2000 Oct 15 


Eyres et al. (2002) 




1998 Jun 29 - 2001 Nov 28 


Skopal et al. (2000, 2002) 




1999 Aug 12 


Eyres et al. (2002) 


UV 


1999 Aug 12 


Eyres et al. (2002) 




1998-2000 Apr 


Eyres et al. (2002) 


X-rays 


- 


- 
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Table A.6. Observation references for the most recent period 
2001-07. 





Observation 


References 




date 




radio 


- 


- 


IR 


1999-2007 


Taranova & Shenavrin (2007) 


optical 


2001 Nov 28 - 2003 Nov 2003 


Skopal et al. (2004) 




2003 Nov 8 - 2007 Jan 15 


Skopal et al. 2007 




2004 Apr 9- 11 


Yoo (2006) 




2004 Oct 2-5 


Yoo (2007) 


UV 


- 


- 


X 


2001 Mar 27 


Galloway & Sokoloski (2004) 




2006 Jan 04/05 - May 28/29 


Mukai et al. (2006) 



